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INTRODUCTION 


Widely different conditions may exist in the ground under buildings erected on 
a concrete slab, depending on the characteristics of the building site, It is there- 
fore difficult to derive universally valid rules relating to the temperatures in the 
ground and to the heat quantities lost through flow into the ground merely on the 
basis of measurements, Mathematical equations have to be set up, which are general - 
ly valid under stipulated boundary conditions, Comparison between the calculated 
values and the results of measurements carried out in experimental houses enables 
rules applicable on a wider basis to be inferred from experimental studies made 
with a few houses only, 

Various theoretical formulae and formulae based on measurements have been 
developed for the purpose of calculating the heat losses through a floor slab, Among 
them, the method rather widely employed in America may be mentioned, in which 
only the losses at the edges are taken into account, The heat quantity given off 
by the slab to the ground ts calculated by the equation 


Q = C(T, - T,) * (2a + 2b), 


where C is a value determined experimentally for various types of structures, 
T; - To ts the temperature difference between outdoor and indoor air, and 2a + b) 
is the circumference of the building [ 1). 

Owing to the entirely different kind of climate and soil types, these experimental 
formulae cannot be used as such in Finnish conditions, 

Billington divides the losses into two components by means of an electric analogy 
method: (1) the edge losses, which consist of the heat passing through the upper 
soil layers into ambient air, and (2) the surface losses, i,e., the heat quantity 
flowing from the central parts of the slab into soil layers at greater depth. Both 
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components are calculated separately in Billington’s method by a simple computing 
procedure, assuming that the isothermal lines are ellipses [ 2], 

Ruckli has devised a calculating method for round slabs, from which the results 
can be transferred to rectangular slabs with the aid of a correction coefficient. 
This method is primarily intended for use in the designing of cold-storage rooms [3 J, 

The methods just mentioned yield mutually different results because they employ 
approximative formulae, 

O, Krischer has presented an approximative method for calculating the heat 
losses going into the ground from buildings erected direct on ground | 4 J; it has 
been modified with a view to more convenient use by W, Weyh| 5 J, Krischer calcul - 
ated the stationary thermal flow from a rectangular building into the ground when 
the constant temperatures Ty and Ty prevail inside and outside and the ground 
temperature is constant at the depth Z = Z,. His solution, derived with the aid of 
series expansions, requires that the heat transfer coefficient (@) be known, How- 
ever, the heat transfer coefficients are dependent on air flow velocity (i.e, on 
the wind) and therefore cannot be considered constant. 

In the present work, the equations of the temperature distribution field and the 
thermal flow conveyed into the ground under stationary conditions are derived 
for a rectangular slab, for a slab shaped as a narrow strip (terrace houses), and for 
a circular foundation slab, These equations are set up for various boundary conditions, 
In order to avoid excessively complicated equations, the said equations in the case 
of non-stationary thermal flow are given only for a narrow strip slab. They 
are given for the following alternatives: (a) The temperature of the slab is a given 
function of time, T(t), (b) the thermal flow from the slab into the ground is a 
given function of time, Qt). 

The temperature of the slab and its variations in terms of time, as well as the 
thermal flow conveyed from the slab to the ground, are obtained from measurements, 
Selecting the temperature function on the lower surface of the slab and on its edges 
in close correspondence with actual conditions, one obtains results conforming with 
measurements, both with regard to the temperature field and the thermal flow, 
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THE TEMPERATURE FIELD PRODUCED IN THE 
GROUND BY A HEATED SLAB 


The temperature in the ground at different depths, 1)(z,t) goes through seasonal 
variations in accordance with the solution of a one-dimensional equation [ 6}. This 
temperature distribution field accountable to the changes of outdoor temperature is 
superimposed upon the temperature field To(x,y,2,t) produced by the slab, When 
the thermal power Qx,y,t) flowing from the slab into the ground is known, the 
temperature field derived from this information can be immediately added to the 
inherent temperature field of the ground, point by point, the resulting, actual 
temperature field being 


T(x,y,2,0) (z,t) + Ty (x, y,2,0 


In the case that the temperature of the slab, T(t), which is assumed to be const- 
ant with respect to the position coordinates, is known as a result of measurements, 
this temperature T(t) already includes the influence of the ground temperature 
field. In calculating the field To(x,y,z,t) produced by the slab, which is super- 
imposed by addition on the ground field, the difference between the actual temper- 
ature of the slab and the ground surface temperature, T(t) - T,(0,t), has to be 
used for the surface temperature of the slab in the boundary conditions, 
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STATIONARY FIELD 


When the ground surface is covered by a snow layer for a length of time and 
the variations of the outdoor temperature are very small, the temperature field in 
the ground can be considered approximately stationary, T, = T,(z), also in the 
upper layers, A slab which does not act as a heat source but obtains its heat supply 
from the room air remains at approximately constant temperature throughout the 
year, and the temperature field under the slab displays rather small variations only. 
Also the field produced by a slab acting as a heat source attains a stationary state 
Ty = To(x,y,2) when the outdoor air temperature ts constant during the time of 
about one month, 


The ultimate field in the stationary state can be described by means of the 
equation 


T = Ty (2) + TolX, 2), 


In the stationary state, both functions appearing in this equation satisfy Laplace's 
equation 


+—+ —— 8 0 (1) 


RECTANGULAR FOUNDATION SLAB 


The temperature outside the slab decreases in accordance with 
a step function, from Ty, to 0, 


The function To9(x,y,z) representing the temperature in the ground, which will 
be written without its sub-index in the following, shall satisfy the differential equa- 
tion (1) and the following boundary conditions: 
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Fig. 1. The choice of coordinates for a 
rectangular foundation slab, 


a) T(x,y,0) * T,, when -a<x<a and “-b<y<b 
b) T(x,y,0) = 0, when -a>x>a and -b> y >b 


9, when -a>x >a and -b? y >b 


d) For Z—poo, 0 


The Fourier double transform of the temperature function 


+e 


1 
T%&,n,2) | T(x, y,2) dxdy (3) 


satisfies the differential equation (1) of the temperature distribution, The object 
function T(x,y,z) is obtained in accordance with the inversion theorem of the 
Fourier transformation [7], 
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Subjecting each term in equation (1) to the Fourier transformation 


+e +e 
2 
| T x+ny) dxdy + €x+ny) dxdy + 
ax” ay? 
+e +e 
| Exeny) dxdy = 0 
82" 


and expanding the terms by partial integration 


+e +e 


+o 
27 
ax2 @x\ax 
+e +o 
( Je dx e ig Ox e NY) dx 


we find upon application of the boundary conditions (2a) and (2b) 


+e +e +e 
o-ig et Exeny dx = 0-ig | -§2 | gy 
+e 


0+0- T(x, y, 2) el Extny) 4, 


+a +e 


ax? 


The second integral of the expression (1), treated in the same manner, gives 
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2 ‘ 
dy? 


10 
&§ 

‘J 
4 

4 

q 

q 

4 

4 
4 

q 

q 


11 


The third integral in (1) is the second derivative of the Fourier transform of 
the temperature function T(x,y,z) with respect to z, 
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By these transformations, equation (1) is converted into a linear, homogeneous 


differential equation of the second order 
") 
( + T* = 0 (5) 
a2" 


The general solution of this equation is | 9 | 


€, s Ve +7 Z + Col et Z (6) 


Only the first term comes into question because the second term increases infinitely 
with increasing Z, When Cy = 0, 


2 
= VE 4 (7) 


We substitute the expression (7) of 1°(€,%,z) in the equation (4), representing 
the temperature distribution, taking into consideration that, according to Euler's 
formula [ 8), 


e "Exeny) cow &x+ny) + i sim &x+ny) 


In this manner we obtain the expression for the temperature field in the form of 
Fourier’s cosine transform 


+e 
1 2 27 
“T(x, y, 2) | Ci cos x+ny)dédn (8) 


when only the real part is considered, for obvious reasons of symmetry. 
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(&,7,9) Ci Qn 


Integration yields further 


T*( €, 9,0) 


The unknown function ce ») in the integrand is determined so that the 
expression of T(x,y,z) satisfies the boundary conditions (2a) and (2b), According 
to equation (6) and (3) and to the boundary condition (2b), 


T(x, y,0) com &x+ny) dxdy = 


+a +b 
&x+ny) dxd 
on COR HY y 


2T,, sin €a sin nb 


En 


8T(x, y, 0) 


dz 


(9) 


(10) 


On the strength of equations (9) and (10), the function T(x,y,z), equation (8), 
describing the temperature distribution acquires the form 


+e +e 
J J sin sin - Ve" + com &x+ny) (11) 


In order that dimensionless quantities might be used, a and b are chosen as 
units of measurement, that is, we write a = @ and nb = #. Bearing symmetry 
in mind at the same time, we now obtain equation (11) in the form 


(12) 


+4) dadp 


Calculating by equation (11), the thermal flow from the slab to the ground, 
per element of surface area, is 


(13) 


and the total thermal flow in the region affected by the slab is 
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When calculated by equation (11), the temperature gradient obtains the value 
oo at the edge of the slab, and consequently the thermal flow will also be infinite 
on the edge. In this case also the total thermal flow from slab to ground becomes 
infinite, as can be seen from its expression (14). 

On physical interpretation, this is found to state merely a self-evident fact: heat 
would have to be supplied to the slab at an infinitely high rate ti order to keep 
its temperature constant up to its very edge when the surroundings of the slab are 
at 0°C, 


This shortcoming of the theory can be corrected by choosing new boundary con- 
ditions, assuming that the temperature of the slab in its marginal parts decreases 
linearly from T,, to zero over a given distance, which depends on the design of 
the slab, This approximation is well consistent with actual conditions encountered 
in practice: owing to the heat release of melting snow and to good insulating 
properties of snow, the ground temperature under the snow cover adjacent to the 
edge of the slab, and also at greater distance, maintains a height close to 0c, 


The temperature in the marginal part of the slab decreases 
linearly from Ty, to 0, 
The new boundary conditions (15) have been chosen on the basis of experience 
gained in connection with actual measurements. 


a) T(x,y,9) = Ty when -a<x<a and -b<y <b 
(X=a) Ty 
b) T(x,y,0) - a<x<a’ 
(y-b) T, 
- "bey <b’ (15) 
d) T(x,y,0) = 0 x >a’ 
e) T(x,y,0) = 0. y>b 


In these formulae, - . mq and - Bes represent the rate of temperature drop in the 
x and y directions, respectively, at the edges of the slab. 

We shall now determine the unknown function C,(€,n) in equation (8) in 
complianee with the new boundary conditions (15) 


+e +e 
1 
T(x, y,0) cos €x+ny)dxdy (16) 
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Fig. 2. Rectangular foundation slab; its temper- 
ature at the edge decreasing linearly from T, 
to zero. 


The Fourier transform of the temperature function on the surface z = 0 is calcul- 
ated as a surface integral, separately for the regions A, B, B’, C and C’. 


sin €a sin nb 


cos §x+ny)dxdy = 4T, 


14 j 

| 

} + 

| 

| 

i 

q 6 

q 

= 

In = Ty | 

| 


15 


7 b 
4 (x-a) Ty 
1, * Ty - co x+ny)dxdy (17) 
a 
a 
(y-b) T, 
b a 
“BY 
In = 0 
= 0 


x (18) 
x com &x+ny)dxdy 


2. 
® afro + f cow § x+ny)dxdy 
a a 
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a 
To ++ Af com Exe - y cox &x+ny) dxdy 
b 
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cos (a’§ n) - cos(a€é (20) 
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b 
a'cos(a’'é +b'm) - acos(aé - bm) sin(a’— +b'm) -sin(ag+bm) 
a a 
sin(a’é -b'n) - sin(ag - bm) (21) 
a 
From equations (20) and (21), we find 
cos (ag +b”) _a cos(ag€ - by) 
a” sin(a’§ -sinag+bn) sina’g -b’n) - 
a’-a +bn)*n (a€ -bn)*n 


The integral I, = I. +1, (equation 15c) fs obtained in a similar manner 
1 2 


c=To + 


(ag (ag +bn)°é 
b? sin(a’é - sin(aé - bn) os 
b’-b (a€ - 


The surface integral for the combination of regions B and C is derived from 


equations (22) and (23) 


cos(agé - cos(a€é - bn) 
- 


( a7 (b’ (a! - a) (sin +b'n) - sin(a€g +bm) 
(a’ - a) (b’ - bag En 


(b’ -b) - b2 (a’ a)) (sin = b’n) sin (a€é (24) 


(a’ - a) (b’ - b) (a€ - 
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For the entire region -a’ < x <a’; -b’ < y <b’, over which the heating 
influence of the slab at ground surface extends, there is according to equations (17) 


sin(a’§ -b'm) - sin (ag -bm) _ 


I I, + +1 = 2T (ae - 


(xt sin(a’é +b'n) - sin(ag +bm) pee 
-b (ab +bn)*én 


As can be seen from Fig. 2, = = a Consequently, the expression in angular 
brackets can be simplified and equation (16) becomes 


ma-a (a€ -by)én (ag +b) En 


The equation of the temperature distribution field is obtained from equations 
(8) and (25) 


+o +e 


T a sin(a’& -b'n) - sin(a€ 
on a'-a f (a& -bn)&n 


(26) 


sin(a’ + b'n) - sin(ag + 


Let us write 


af =@ 
bn = 8B 


b 
and hence ag 
a 


by which the solution can be presented in the following form, which is more 
appropriate for numerical calculations 


+ote 
2Toa sin a’ 
T(x,y,2) xia a_(@-B) - sin( _ 


(a -B)aB 
sin 2’ a Ji (22) cos * cos dadB (27) 


This equation can be employed to calculate the temperature field in any point (x, 
y,Z) under a rectangular floor slab with the sides a and b, produced by this slab 
in the stationary state. 
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In principle, the calculation can be carried out by giving B, at first, a fixed 
value B, and performing numerical integration with respect to @, using tables 
[10]. All values of # are gone through in this manner, the result being that 

F (a@B,) d@ is obtained as a function of B. Subsequently, graphical integr- 
ation with respect to B gives the temperature in the desired point (x,y,z). 

Even with the aid of tables, the calculation involves a very great amount of 
work and it is virtually necessary to resort to a computer. (The IBM 710 computer 
was used in the present instance.) 

We are now going to show that, for b-»oo, the temperature function (27) ap- 
proaches as its limit the solution (45), p. 25, obtained for a slab having the shape 
of a narrow strip. 


Let I, denote the inner integral with B as its integration variable 


(a-B)aBp (a+B)aB 
Xa + cos 
lim cos + xB) cos 
when b —> oo 
2 
a Zz 
lim e “ae + (58) ze 


when b °° 


As its limiting value, the integral I, obtains the form ly = I + Ip + I, + I 


sin (@-B) _ sin (@ -B) sin dp 
ha (a-B)aB (a+p)ap 


Each integral in the expression of I, is calculated separately, the results being 


+o 
sin (@ % cos @ 
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a 
(@+f) n 


(a+B)aB 


sin (@ +8) cosa 


Summation of the four integral renders 
and hence further 
+% cos @ cos e cosa 
a 
The limit of the temperature function (27) for b—» oe will thus be 


+e a’ 
2a Ty cos — - cos @ 
m(a'-a) J a 


Zz 
T (x,y) = - e cos da, (28) 


which is the same as equation (45), p. 25 (a = b). 


The thermal flow from a rectangular slab to the ground. In 
order to maintain in the ground under the slab a stationary temperature field in 
accordance with equation (26), the heat quantity 

aT 
dQ=-4 dxdy (29) 
per surface area element has to be conveyed from the slab to the ground. The 
temperature gradient on the lower surface of the slab is obtained from equation 
(26) by derivation with respect to z 
+e +e 
aT (x,y, 0) a -b'n) - sin (ag 


az (a€ -bn) En 


sin (a'§ +b'n) - sin (a§ +bn) 
+ bnyén + cos &x+ cos ny (39) 


{ 
i 
} 
4 
: 
i 
j 
| 
| 
| 
i 


20 


Substituting the expression of this gradient in equation (29) and carrying out the 
necessary integration, we find the tota) thermal flow 


+e+e »b 
0 0 -a 
AT, -b'n) - sin (ag -bn) 
a sin (a -b 7) - sin (ag -by) 
| (ag -by)&n 


sin (a'€ +b’) -sin (a sin Mb 2 
déd 31 


The numerical values involved in equation (31) are calculated in like manner as 
those in equation (27). 


Using the boundary conditions (15), physically correct solutions are obtained for 
the temperature field as well as the thermal flow. This will also be shown with 


the aid of numerical calculations, applying the equations derived under b) for a slab 
having the shape of a narrow strip. 


FOUNDATION SLAB HAVING THE SHAPE OF A NARROW STRIP 
Ground surface temperature outside the slab is 0°C. 


If a >b (elongated buildings, terrace houses, etc.), the temperature field in 
the middle of the slab can be calculated by assuming the slab to be a rectangle 
of infinite length and of width a. 


The temperature function T(x,z) shall satisfy Laplace’s equation 


a2T 
(32) 


and the boundary conditions 


a) T(x,0) = when -a<x < ta 
b) T(x,0) = 0 7 -arx > ta (33) 
c) T=0 for z = co 
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Fig. 3. Foundation slab having the shape of a 


narrow strip. Temperature of the slab T,; tem- 
perature at ground surface outside the slab 0°C. 


The Fourier transform 


+o 
= ax (34) 
satisfies equation (32). 
According to Fourier’s inversion theorem, 
+e 
T (x,2) = > ge (35) 


Equation (32) can be solved by multiplying it by e 1ez Partial integration con- 
verts equation (32) into 


Qe 
= T 
(36) 


‘ Selecting the finite branch of the exponential constituting the solution of equation 
) (35), we find according to boundary condition (33c) 


T% = c(€) 
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22 
By virtue of boundary conditions (33a) and (33b) and employing Euler’s formula 
= cos +i sin &x (38) 


we have, when the real part of the function is chosen, 


8,0) =c(8) V2 con bx dx = 


According to equation (35), 


T(x, 2) ame (39) 
0 


which is a function satisfying the boundary conditions (33). 
Let us now introduce dimensionless variables by employing a as a unit of measure, 
and writing a = @. Equation (37) then becomes 


Zz 
T(x, 2) cosa e* da arc (40) 


0 


The graphical representation of equation (40) is shown in Fig. 4. 

The same result is obtained from the theory of the rectangular slab if one of 
its sides is allowed to increase infinitely. This is easily seen by determining the 
limit value which is approached when b—»oo in equation (43). 

In deriving equation (39), it was assumed that the temperature decreases from 
T, to zero in the manner of a step function at the edge of the slab, This is not 
consistent with actual conditions, the decrease of temperature in the lagging 
on the edge of the slab actually following a certain continuous function. If 


we approximate this function by means of a linear function, we find new bound- 
ary conditions: 


a) T(x,0) = Ty when -a<x<a 

X-a 
b) 27, a<x<b (41) 
-*- £06 x x >b 


d) for 
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Fig. 4. Course of the isothermal lines and thermal flow from slab to ground 
under a foundation slab having the shape of a narrow strip according to 


equation (40) when T, = 30°C. 


We shall determine the unknown function C(§) in equation (37) in conformity 


€,0) =C(€) cos &x dx 


with the new boundary conditions (41) 


b Ex dx =F 


0 


(42) 


Integration, considering the boundary conditions (41a), (41b) and (41c), renders: 
- in the interval 0< x <a 


(b - a) sin €a 
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Fig. 5. Foundation slab having the shape of a 
narrow strip; its temperature at the edge de- 
creasing linearly from T, to zero. 


- in the interval a < x < b 


a 


- and for x > b 
Il, = 0 


T, cos &b - cos éa 
(a - b) 


The Fourier transform for the boundary conditions on the surface will thus be 


2 To cos - cos €a 
T% €,0) = C(€) (43) 


According to equations (35), (37), (38) and (43), the temperature function satis- 
fying boundary conditions (41) will be 


2 Ty cos §b 
T(x, 2) * 
0 
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Finally, introduction of dimensionless variables by writing €a = @ converts equation 
(44) into a suitable form for numerical calculations 


cos—@ - cos @ a 


The graphical representation of (45) is shown in Fig. 6. 


The thermal flow into the ground.The thermal flow production under 
the slab a temperature field in accordance with equation (43) is obtained from the 
equation 


@T(x, 0) 


dQ=-A 


(46) 


The derivative of the temperature function with respect to z at z = 0 gives the 
gradient at the lower surface of the slab 


OT(x,0) 2T, {= - cos 


az m(a - d) 


The total thermal flow is 


m(a-b) 


4AT, f (cos &b - cos sin &a 
m(a - b) he 


dg (47) 


The temperature field and thermal flow for a slab having the shape of a narrow 
strip, calculated by equations (45) and (47), mm are shown in Fig. 6. They are 
well consistent with the results of measurements in homogeneous ground when the 
temperature field and thermal flow inherent in the ground are superimposed on 
those shown by Fig. 6. 
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Fig. 6. Course of the isothermal lines and thermal flow from slab to ground 
under a foundation slab having the shape of a narrow strip according to 
equation (45) when T, = 30°C. The downward heat flow is positive. 
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The integral 


Ba cos &x dx dé dé (48) 


representing the total thermal flow calculated from the temperature function (39) 
derived on the basis of the boundary conditions (33) is not finite. Integration with 


respect to x from 0 to x gives a formula for the thermal flow in a zone from 
the centerline to a line at the distance x 


= 


This equation can only be used in the central part of the slab. 
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The temperature function T(x,z) defined by equation (44), p. 24 approaches 


as its limit the temperature function in equation (39), p. 22 when b approaches a. 
In order to show this, we write 


b-a 


lim oe = 2 lim sin 3 a 


cos &b - cos &a &(b + a) ‘ein sin 


b—a 
= 2 sin = sin 


by which operation the equation (44) is transformed into 


which is the same as equation (39). The consideration shows that the calculations 
have been correctly performed. 


CIRCULAR FOUNDATION SLAB 


In the following, the stationary-state temperature field and thermal flow under 


a circular foundation slab shall be considered, choosing the boundary conditions in 
three different ways. 


The surroundings of the slab are completely insulated 


We introduce the cylindrical coordinates r, g, z. The temperature is inde- 
pendent of @ and Laplace's expression has the form 


1 aT 


This is a Besselian differential equation of the order zero. 
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The temperature function T(r,z) shall satisfy the equation (50) and the following 
boundary conditions: 


a) T(r,0) = Ty when r<R 
b) = 0 r>R 
c) 20 r= 0 
d) T(r,z) 0 oo 
oo 
The Hankel transform of the order zero of T(r, z) 
T%&,z) = r T(r,z) Jo (€ndr, (52) 


satisfies equation (50). 
We subject equation (50) to Hankel’s transformation: 


Oz 


2 
T(r,z) 1 @T(r,z2) 
Be r 


aT 
= — J Ende 


The Besselian differential equation (50) is converted into a linear, homogeneous 
differential equation of the second order, 


2 


az 


The Hankel transform of the temperature function thus satisfies the differential 
equation (53). 


The general solution of equation (53) is 
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In compliance with boundary condition (5lc), we choose the finite branch of the 
function, which implies that B(€) = 0 


T% = (gy (55) 


The unknown function A(€) in this equation will be determined by the boundary 
conditions. From equations (52) and (55), we obtain 


According to Hankel’s inversion theorem, we have further 


T(r,z) EA Ee Ende (56) 


Differentiating equation (56) with respect to z, we find 


Applying the boundary conditions (51a) and (51b) to equations (56) and (57), we 
obtain the dual equations 


T(r,0) EAE =T, 
(58) 
I, 2 
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Choosing the radius R of the circular slab as a unit of measure and writing 


t u 
—=@, §=— 
0o<e<il 
we have 
FACE) I = Jz u (2) = Ty 
3 If = u2 = 0 
5 


Let be u? a(#) = f(u); then u A(z) =u f(u), and the dual integral equa- 
tions (57) are converted into 


. 
J fu) = R27, 0<e<1 


(59) 
f(u) = 0 @>1 


The function f(u) is determined with the aid of these two equations, 
The equations just obtained constitute a particular case of the more generally 
valid dual integral equations 


u® f(u) J, (ug)du = ge) 0<e<1 
(60) 
f(u) J, (ug)du = 0 
where g(@) is a known function and f(u) is being sought. 
The dual integral equations were first solved by Titmarsch [11] with the aid 


of Mellin’s transformation. His solution is only valid if @ >0. The equations 
have also been solved by I.W. Busbridge, his solution being valid when @ > -2 and 
“9- 1<(« -4)< ¥ +1, provided that g(u) is integrable in the interval 0 to 1 [12]. 
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Fig. 7. Circular foundation slab. Constant tem- 
perature of the slab = T,,; surroundings of the 
slab completely insulated, 0, 


Busbridge’s solution is 


f(@) = 


Jy ot. g(@) for? (1-2)? gu)du + 


Oo 


1 1 

sa 
for (a - v4" dy J cew Jy 41 + (@u)du 


In the case that @ = -1 and » = 0, and g(u) = R°T 


o» the unknown function 
in the equations (59) is 


f(@) (@) fo au + 


( eu?” (@u)du (61) 
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We substitute in equation (61) 


r(z) 
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1 1 
2 
Jy = (@) 11 Vag cove 
2 


2 


2 /sin@ 
Ign (@) = (“5° = cos e) 
by which this equation is converted to 
2 Re (62) 
R°T V2 2 
The unknown function in the equations (59) is thus 
QR°T 
f(u) = — sin u (63) 


and the unknown function A(€) is obtained by re-introducing the initial variables 
u = RE 


The final form of the temperatuse function (56) will then be 


2T R 
re) 


On the boundary surface z = 0, the integral occurring in the solution becomes 
% for = R and arc sin for r > R (13). 

The function T(r,0) equals for r= R and T(t,0) = 21% 
arc sin R for r > R. When t—moo, T(r,0) approaches the limiting value 0, The 
stipulated boundary conditions (51a) and (51c) are thus satisfied. 


There is 
, RE =0 when r > R 
az n 


and equation (65) thus satisfies also the boundary condition (51b). In the solution 
outlined here, the unknown function was sought by applying a systematic procedure. 
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Fig. 8. Course of the isothermal lines and thermal flow from slab to ground 
under a circular foundation slab according to equation (66) when T, = 20°C. 
The downward heat flow is positive. 


It can also be found by trial [14]. For the purpose of numerical calculations, 
we write R& = @, i.e., we choose the radius of the circle, R, as a unit of 
length. The integral will then be 


Z 
2T ina 
T(r,z) = —2 e ' z>0 (66) 
a R 


The temperature field defined by equation (66) is shown in Fig. 8. When the ground 
surface around a circular slab is completely insulated (boundary condition 2), the 
isothermal lines are perpendicular to the ground surface outside the slab. This 
condition is approximated in actual practice when the ground surface is covered 
by a thick layer of snow. 
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The thermal flow into the ground, In order to maintain a stationary 
temperature field in accordance with equation (65) under a circular floor slab, the 
thermal flow from the slab into the ground which passes through an annular element 
of width dr has to be 


OT(r, 


dQ -24 ——— dr (67) 
According to equation (65), 
0 


The integral obtains the following values [ 13; p. 405 ] 


1 
sin RE dE = when OC 


2 (0 r>R 
= or 0 Rer 


Equation (67) is thus converted into 


dQ = (68) 


The total heat quantity flowing from the slab into the ground is obtained by inte- 
gration: 


R 
Q=+4T,A = 4T AR (69) 


It can be seen from this equation that the total thermal flow from the slab into 
the ground is directly proportional to the thermal conductivity of the soil, A, and 
to the radius of the slab, R. As the floor area of a building increases with the 
square of this radius, it can be inferred that large buildings are more advantageous 
than small buildings from the viewpoint of heat losses. 

Equation (68) reveals that the losses increase rapidly towards the edge (when 
r-*R) and become infinite on the boundary line, while the total thermal flow 
itself remains finite (Fig. 8). 
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Ground surface temperature outside the slab is 0°C, 


We proceed to solve the Besselian equation (50) under the boundary conditions 


a) T(r,0) T, when 
b) T(r, 0) = 0 r>rR (70) 


As in the preceding instance, Hankel’s transformation and Hankel’ s inversion theorem 
can be employed to obtain the temperature function 


According to the boundary conditions (70), equation (71) renders the dual integral 
equations 


T(r,0) = = T, O<r<R 

(72) 

T(r,0) = = 0 r>R 

As a new unit, we introduce the radius of the circular slab, writing 
and when € increases from 0 to oo 


0<e<i1 then u increases from 0 to —® oo 
= gu ua(?)- f(u) 
The dual integral equations (72) now become 
flu) @u)du = 0<e<1 


(73) 
Jf e@uddu = 0 


These equations represent another particular case of the dual integral equations 
(60) (for = 0, ¥ = 0 and g(u) = RT,). 
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In this case, too, the solution of equations (72) is obtained from Busbridge’ s 
solution by substituting @ = 0, ¥ = 0 and g@) « PT, for g(@) is integrable 
in the interval 0 to 1 and regular for y=®ee, Moreover, the conditions @ > - 2 
and -»-1<{@ are satisfied, 


1 1 1 
fe) = du + v dv f eu)” J,( @uddu = 


1 
1 


We introduce a new variable, writing 
du 4 dt 


=t, u® 


and carry out the integration 


= @Jole@) 


1 

1 
(gu) = 7 


@ 
J,(odt = Jo @) 

Hence 


1 1 
f(@) = RT, [ este) + ese) | = 5R°T, e| 
1 2 


fig) = RT, (14) 


Re -introducing the initial variables 


f(u) = J(u) 
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J,(R€) 
Ag) = 


we find the final form of the temperature function (71) 


= 
T(r,2) = RT, Joft €) dg (15) 


At ground surface (for z = 0) this temperature function has the values 


O whenr>R 


T(t,0) RT, J(RE) dé = "pf 
To 


With a view to easier numerical calculations, we choose the radius of the slab 
for the unit of measurement and write RE =a@; Rdg = da; We have 
then 


T(t,2) Ty toler) e da (76) 


The isothermal lines of the temperature field consistent with the boundary con- 
ditions (70), calculated by equation (76), have been shown in Fig. 9. 
The temperature gradient calculated by equation (71) for z = 0 


A T(r, 0) 


and the total thermal flow 


Q = 2n ART, fi (RE)dE 


become infinite. No physically correct result is thus obtained under the boundary 
conditions (70). 
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Fig. 9. Course of the isothermal lines and thermal flow from slab to ground 
under a circular foundation slab according to equation (76) when T, = 20°C, 
The downward heat flow is positive. 


Constant thermal flow from the slab, and completely insu- 
lated surroundings 
Let a constant thermal flow Q, be conveyed into the ground through the area 
under a circular slab, assuming that the thermal flow at ground surface outside the 
slab is zero. 


Owing to symmetry of the temperature field, Laplace's expression in polar 
coordinates is not dependent on @, and it has the form 


e2r 10T 
am * or * 


which is a Besselian differential equation of the order zero. The temperature function 
T(r,z) shall satisfy equation (77) and the boundary conditions 
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Fig. 10, Circular foundation slab; the thermal 

flow Q, entering the ground from the slab 

being known. Completely insulated surroundings 
of the slab, 


oo 

co 

Hankel’s transform of the order zero of the temperature function T(r,z) for the 
derivatives, 


d) T(r, 2)» 


€,z) = t T(t,2) z>0 (79) 
satisfies equation (77). 
We subject equation (77) to Hankel’s transformation 


&, z) 8“T(r, z) z) z) 


and expand the integral by partial integration. Taking into account the boundary 
conditions (78c) and (78d), the Besselian equation (77) is then found to be converted 
into a linear, homogeneous equation of the second order 


,z) 


= £2768 
(80) 


which is satisfied by Hankel’s transform (79) of the temperature function. 
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The general solution of equation (80) is 


T% ac eye” + wg) (81) 


By virtue of boundary condition (78¢), B ~ 0 and the solution has the form 


We shall then determine the unknown function A(€) in equation (82) on the 
basis of boundary conditions (78a) and (78b). 


According to Hankel’s inversion theorem, we obtain from equation (79) the tem- 
perature function 


T(r,z) &,z) If Ende (83) 
0 


Equation (82) can be used to convert the temperature function (83) into 


T(t,z) fe yy (84) 
and we find by derivation 


and according to boundary conditions (78a) and (78b) in the plane z = 0 


fe ACE) If 7 % O<r<R 
ACE) If = 0 
az 


The unknown function A(é) is determined from the dual integral equations 


fe ACE) Jf Edge = % Oo<r<R 

(85) 
= 0 r>R 
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We write eg? A(8) = GE), by which the equations (85) are converted into 


- 
fxn O<r eR 
fae) If = 0 Rer 
0 


This is a particular case of the dual equations (60), p. 30, for a= 0 and v ® 0 


r u = 
= gu 
f(t) Q 
folk = R= 


o(i:) Jo @uddu = 0 


| We write o(it) = flu); 


= 
(86) 
f(u) Jof@u)du = 0 
| 
| The solution of the dual equations is found to be (cf. p. 36) 
Q 
f(@) = J(@) 


f(u) = J4(u) 
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1 RQ, 
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Fig. 11. Course of the isothermal lines under a circular foundation slab, 

and surface temperature as a function of the radius, according to equa- 

tion (87), when the thermal flow entering the ground from the slab is 
Qo = 6,54 kcal/m?h. 


~ 


Substituting the expression for A(€) in equation (84), the temperature function 
is found to be (Fig. 11) 


which satisfies the stipulated boundary conditions. (Cf. Weber-Schaftheitl's integral 
[13; p. 398). 
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NON-STATIONARY FIELD 


SLAB HAVING THE SHAPE OF A NARROW STRIP; ITS SURFACE TEMPERATURE 
VARYING AS A FUNCTION OF TIME 


In the cases dealt with in the preceding chapter, the temperature of the heated 
slab was assumed to be constant and the temperature field in a stationary state 
was considered. However, the temperature of a foundation slab varies during the 
heating period according to the heat requirements of the building, depending on 
weather conditions, and it follows the temperature changes of ambient air. The 
thermal flow from the slab into the ground, Q (in keal per m? and hour) is also 
a function of time, which is known on the basis of measurements 


Outside the slab, along the wall of the building, the snow that has fallen during 
the winter, together with snow coming down from the roof, forms a thick layer 
of high thermal insulating capacity. The temperature at ground surface is therefore 
close to 0°C. In the mathematical treatment it will be assumed to equal zero. 

The temperature field produced by the heated slab in the ground is a periodic 
function of time T(y,z,t), upon which the periodic temperature field of the ground 
itself, T(z,t), is superimposed, 


The function To(y,z,t) representing the temperature field produced by the slab 
satisfies the differential equation 


(89) 
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Fig. 12, Foundation slab having the shape of a 
narrow strip; its temperature being a given func- 
tion of time, T(t). 


if the ground is assumed to consist of a homogeneous, isotropic substance having 
constant temperature conductivity a. We solve this equation under the boundary 


conditions 
T(y,z,0) = 0 (90) 


and 


a) T(y,0,t) = T(y,t) = T(t) when (-b < y < b) 


« 0 /y/> b 
b) Fy 10.2.9 = 0 (91) 


c) T(y,z,t)—>0 when y and z —> oo 


In order to find the solution of equation (89), it is converted into a linear, homo- 
geneous differential equation by means of Fourier’s sine and cosine transformations 


2 
me 


J t.20 cos ny sin Oz dy dz (92) 


To this purpose, both members of equation (89) are multiplied by 4 cos ny sin Oz 
and integrated 
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ny sin Oz dydz = + ny sin Oz dydz 
at 822 

00 


ee 


t 


at 


a4 
+b i 
i 
a 
| 
| 
| 
| 
i 
rey | 
| 
| 
aT 
4 
4 


45 


2a 

abe cos ny sin Oz dy dz = 
0 


a 2a 2 
CT(y,0, nydy - a ff I(y,z,theos ny sin Oz dydz 


Second integral: Taking into account boundary conditions (91b) and (91c), we find 


Equation (89) is thus converted into a linear, homogeneous differential equation of 
the first order 


cos my sin Cz dydz = - a 0,0 


+ a(n’ 0,0 = 0 (93) 


Its solution can be obtained, e.g., by the method of variable constants | 9; p. 49 ] 
t t t 
2472 
n 


(94) 


t 
at(n“+ a 
According to boundary equation (90), 
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T% 9, 0,0) 


cos ny cos dydz = 0 


The solution (93) of the linear equation has therefore the form 
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The solution fo the initial equation (89) is obtained from equations (92) and (95) 
by means of Fourier’s inversion theorem 


eat (96) 
2 
The integral involved in this solution is Fourier’s sine transform of the function 
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32 
4(a(t-t) 


The temperature function (96) - thus be written 
2 
T(y,2,) cos my dn dt (97) 
We shall determine the function Te(n, t) by the boundary conditions (91a) 


Tg (n,t) T(y, cos ny dy 


According to boundary condition (91a), the function T(y,t) obtains the value T(t) 


when -b < y < +b, and it is zero when /y/ > b. Consequently, we find by inte- 
gration 


Inserting the expression (98) for TZ in equation (97), the temperature function is 
found to be 
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The function G,(,y) * “a is Fourier’s cosine transform of the function 


-/2 when -0 < y < 4b 
} 
when /y/ = b 


2 
and Gon, 4s Fourier’s cosine transform of 
4a(t- 


Vi 


The integral Gy(m,y) Gol, y)cos nydn is Fourier’s cosine transform of the 
product of two Fourier transforms. By convolution [7; p. 44), we obtain 


+v)2 

Substituting the integral in equation (99), the solution is found to be 

_. 2 2 
lt @ 4a(t- 7) (y+v)" 
a ze 4a(t-T) 4a(t-T) 
T 
(y,2,0) = 4a e +e (t)dvdt (100) 
00 

| We shall further modify this solution for more convenient numerical calculation, 


writing 
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which converts the inner integral into 
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Fig. 13. Surface temperature of the slab, and thermal flow that has entered 
the ground as a function of time measured in one of the experimental houses. 
Starting point: lst of September. 


The integral from b to co vanishes 


-b y+b 
4 4a(t-T) 
-V4a(t-t) e dp + V4a(t-t) e Y dy = i 
y 
4a(t-T) 


- y+b y-b 
yz Vv 4a(t Tt) erf 


Substituting the resulting expression of the integral in equation (100), we have 
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(101) 
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Fig. 14. Course of the isothermal lines under a foundation slab having 

the shape of a narrow strip, and its surface temperature at the time 

t = 200 days when its surface temperature varies as a function of time 
in conformity with Fig. 13; a = 0,0032 m?/h. 


Equation (101) represents the temperature, due to the effect from the slab, 
in the point (y,z) at the time t, when z > 0. When T(t) is known 
(Fig. 13), the temperature field is obtained by means of equation (14) 
(Fig. 14). 

The expression of the function reveals immediately that it satisfies the boundary 
conditions (90) and (91a), (91b) and (91c). It is also easy to compute with the 
aid of tables and graphical integration. 

By the limit process t— eo, equation (39), p. 22, describing the temperature 
distribution in the stationary state, is obtained from equation (101). In order to 
show this, we shall return to equation (100) and integrate it with respect to t, 
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considering, at first, that lim T(t) = Tg when too. From equation (100) 
follows 


T 
T(y,z,0) e 4a(t-T) +e 4a(t-T) dvdt 
4 4a(t- 


Tt) 
(y-w? + (yn)? 
+ (y-v)2 + (yw 
2” + (y-v)? 
v) 22 + (y+v)" +v)" 
Zz 4at 
n + (y-v) z + (y+v) 


The limit value of the temperature function for t-—»a@_ is found to be 


Zz 
= 


The integrand is Fourier’s cosine transform of the function 


2 
cos (yn) e 


and the temperature function is 


T 2 2 
T(y,0 [co (y n) cos (vn) e andy 
0 


We integrate with respect to v 


2 


oO 


cos yn e dndv 


This is the same equation as (39) on page 22, if we replace y by x and by €. 
The thermal flow into the ground has been calculated from the isotherm field 
(Fig. 14). 
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SLAB HAVING THE SHAPE OF A NARROW STRIP; RELEASING INTO THE GROUND 
A THERMAL FLOW VARYING AS A FUNCTION OF TIME 


In the following, we shall consider the case that the thermal flow conveyed 
from the slab into the ground, Q (in kcal per m*” and hour), is a function of 
time, known from measurements. 


oT 
= (103) 


As a result of the thermal insulation provided by the snow cover, the thermal 
flow from the ground surface into the air is small; it shall be assumed to equal 
zero in the mathematical treatment of the problem. 


The function Toly,2,0) representing the temperature field produced by the slab 
satisfies the differential equation 


a ay? + (104) 


We solve this equation under the boundary conditions 


T(y,z,0) = 0 (105) 


and 


= - 1 - 
a) T(y,0,t) = Qy,t) At) when -b<y<b 
= 0 bey<-b 


b) T(0,z,t) = 0 (106) 


c) T(y,z,t)—» 0 when y and oo 


In order to find the solution, we convert this equation into a linear, homogeneous 
differential equation with the aid of Fourier’s double transformation 


[ cos ny cos Cz dy dz (107) 
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Fig. 15. Foundation slab having the shape of a 
narrow strip; the thermal flow entering the ground 
from the slab being a given function of time. 


For this purpose, both members of equation (104) are multiplied by & cos ny cos Oz 
and integrated, considering the boundary conditions 
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Equation (104) thus obtains the form of a linear, homogeneous differential equation 
of the first order 


tun? + t% = 0 (108) 


Its solution is 


T*(, 0,0) + 


t t H 

v2 


t 
-at 2 
According to boundary condition (105), 
ee 
0,0) cos my cos Oz dydz = 0 
0 


The solution (109) of the linear equation becomes therefore 


t 
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The solution of the initial equation (104) is obtained from equations (107) and (110) 
according to Fourier’s inversion theorem 


eet 2 (111) 
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The integral involved in this solution is Fourier’s cosine transform of the function 
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The temperature function can thus be written 


22 


(112) 
2 

V2 e 

Ty,2,0 cos ny dn dt 


We determine the function Q*7,t) by the boundary conditions (106a) 


2 
“Vz (y,t) cos ny dy 


The function Qy,t) obtains the value-5 Qt) when 0 < y < +b and the value 
0 when /y/ > b. This function has the Fourier transform 


The temperature function (112) can thus be written 


z2 (113) 
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The function g,(%,y) au is Fourier’s cosine transform of the function 
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The integral ycos nydn is Fourier’s cosine transform of the 
product of two cosine transforms, By convolution, we obtain 


» Van (y-v)? (y+)? 

4alt-t) da(t-t) 


Substituting the integral in equation (113), the solution is found to be 
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T(y,z,0 = + QA t)dvdt 


We shall further modify this solution for more convenient calculation, writing 
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Substituting the resulting expression of the integral in equation (114), we have 


ssss 


T (yo.t) 


9.73 
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Fig. 16. Course of the isothermal lines under a foundation slab having 
the shape of a narrow strip, and its surface temperature as a function 
of y at the time t = 200 days when the thermal flow varies as a 
function of time in conformity with Fig. 13; a = 0,0032 m@A, 
A = 1,00 kcal/n h°C. 
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T(y,z,t) = Wi 
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Q(t) is known from measurements (Fig. 12). The temperature field produced by 


the slab is obtained from equation (115) (Figs. 16 and 17). 
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Fig. 17. Course of the isothermal lines under a foundation slab having 
the shape of a narrow strip, and its surface temperature as a function of 
y at the time t = 200 days when the thermal flow varies as a function of 
time in conformity with Fig.13; a = 0,00022 m@h, A = 1,43 keal/m h°C, 
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SUMMARY 


The mathematical equations of the temperature distribution field produced by 
a warm floor slab in the underlying homogeneous, isotropic ground have been 
solved and the heat quantities flowing from the slab to the ground, in the station- 
ary state, have been calculated for slabs of rectangular form, of the shape of a 
narrow strip, and of circular shape. In the solutions those boundary conditions have 
been found which are most appropriate for use in cases occurring in actual practice. 

Non-stationary temperature fields for which the surface temperature of the slab 
or the thermal flow from the slab to the ground is known as a function of time 
are presented for a slab having the shape of a narrow strip. The equations can also 
be employed to calculate the field, varying with time, which exists in the vertical 
section in the middle of an elongated, rectangular slab. Fourier’s and Hankel’s 
transformations have been employed in the solutions and the unknown functions in- 
volved have been found by means of a systematic procedure, e.g., as solutions 
of dual integral equations. 


| 
| 
| 
ee 
| 
| 


REFERENCES 
Heating, ventilating, air conditioning guide: 1956. Baltimore 1956. 


Billington, N.S.: Electrical models for heating problems. Heating and 
ventilating engineer 24 (1950):282 p. 245 « 247. 


Ruckli, Robert: Der Frost im Baugrund. Wien 1950. 


Krischer, O.: Die W&armeaufnahme der Grundfldchen nicht unterhellerter 
Raume. Gesundheitsingenieur 57 (1934) S, 513. 


Weyh, W.: Die Berechnung des W4rmeaustausches von Bodenfldchen ge- 
heizter oder gekilhlter Raume. W4rme- und Kadltetechnik 38 (1936):11. 


Vuorelainen, O.: Thermal conditions in the ground from the viewpoint 
of foundation work, heating and plumbing installations and draining. Helsinki 


1960. [ Acta Polytechnica Scandinavica, Civil Engineering and building con- 
struction series 7 


Sneddon, lan N.: Fourier transforms p. 105. New York 1951. 
Lindeléf, Ernst: Johdatus funktioteoriaan s. 63. Helsinki 1936. 


Lindeléf, Ernst: Differentiaali ja integraalilasku ja sen sovellutukset III, 
s. 132. Helsinki 1935. 


10. Jahnke-Emde: Tafeln héher Funktionen. Leipzig 1948. 


| 
a 
2. 
| : 
6. 
1 
| | 
8. 
3 


60 
11. Titmarsch: Theory of Fourier integrals p. 335 ~ 337. Oxford 1937. 


12. Busbridge: Dual integral equations, London 1938. [ Proceedings of the London 
Mathematical Society 44, p. 115). 


13. Watson, G.N.: A treative on the theory of Bessel functions, 2. ed. p. 405. 
Cambridge 1958. 


14. Carslaw, H.S, & Jaeger, J.C.: Conduction of heat in solids. p. 188, 
Oxford 1947. 


- 
i 
| 
| 
| 
| 
| 
ay 
j 
| 
| 
‘ 
| 
| 
. 


| 
| 

| 
} 
| 
| 
agt 
| 
¢ 
| 
\ | 
a 
} 

{ 
| 


ae 

ay 

& 

ag 


4 

4 

a 

ives 


THE LAST VOLUMES OF ACTA POLYTECHNICA 
Civil Engineering and Building Construction Series 
(The predecessor of Acta Polytechnica Scandinavica) 


Volume 4 


Nr 1 Lanoz Hansen, P: The Plastic Theory of Curved Beams with Compressive Axial Forces. Acta P 215 
(1957), 32 pp, Sw. Kr 5: 00 UDC 624.072.9:539.414 
Nr 2 Isaksson, A: Creep Rates of Excentrically Loaded Test Pieces. Acta P 219 (1957), 33 pp,.S™ Kr 4: 50 
620.172.251.4 
Nr 3 Iporn, G M: Concrete Deterioration of a Foundation. Acta P 221 (1957), 48 pp, Sw. Kr 8: 00 
UDC 620.163.24:693.5 
Nr 4 Orpino, F Bove: Die Bestandigheit der Basislatien,. Acta P 225 (1957), 16 pp, Sw. Kr Zi vy P 
526.495 
Nr 5 Sunpstrrém, E: Creep Buchling of Cylindrical Shells, Acta P 230 (1957), 44 pp. Sw. Kr 5: 00 
IDC 
Nr 6 Kupsx-Jorornsen, B: Calculation of Time-Concentration Curves in a Stationary, Laminar Liquid 
Flow Through a Circular-Cylindrical Tube. Acta P 233 (1957), 18 pp, Sw. Kr 5: 00 
UDC 
Ne 7 Enozziunp, F: On the Theory of Multiple-Well Systems. Acta P 234 (1957), 11 pp, Ov. Kr 5: 00 
551-495-34 
Nr 8 Baincn Hansen, J: Calculation of Settlements by Means of Pore Pressure Coefficients. Acta P 235 
(1957), 14 pp, Sw. Kr 5: 00 UDC 624.131.5326 
Ne 9 Daut, N J: On Water Supply from Wells. Acta P 236 (1957), 17 pp, Sw. Kr 5: 00 
UDC 551.49:628.112 
Nr 10 Lunporen, H: Dimensional Analysis in Soil Mechanics. Acta P 237 (1957), 32 pp, Sw. Ks ¢ 6: 00 
UDC 624.151 


ACTA POLYTECHNICA SCANDINAVICA 
Civil Engineering and Building Construction Series 
Ci 1 Bretrino, A E: Stable Channels. (Acta P 245/1958), 130 pp, Sw. Kr 7: 00 UDC 626.01:67.1 


Ci 2 Osterman, J: Notes on the Shearing Resistance of Soft Clays. (Acta P 263/1959) 24 pp, Sw. Kr. 7.00 
UDC 624.131.222 


Ci 3 Eme, Owe anp JOHANNgSEN, Ivan J: Measurement of Strut Loads in the Excavation for Oslo Technical 
School. (Acta P 266/1960) 14 pp. Sw. Kr. 7.00 UDC 624.131.532 : 624.1344 


Ci 4 JOHANNESEN, Test Section and Installation of Test Oslo Subway and Kyetr: 
An Earth Pressure Cell for use on Sheet Piles (Acta P 267/1960) 16 pp, Sw. Kr. 7.00 
UDC 624.131.3686: 624.1344 
ct 5S Kyaansut, Byorn: Test results, Oslo Subway (Acta P. 268/1960) 11 pp. Sw. Kr. 7.00 
UDC 624.131.532:624.134.4 


Ci 6 Kyrxtunp, Haratp: Uber die Einschatzung von Biegespannungen in gekriimmten Balken. (Acta P, 274) 
1960) 56 pp. + fig. Sw. Kr. 7:00 UDC 624.072.7 
Ci 7 Vuorerainen, O: Thermal Conditions in the Ground from the viewpoint of Foundation Work, Heating 
and Plumbing Installations and Draining. (Acta P. 277/1960) 40 pp. Sw. Kr. 7: 00 
UDC 624,131.43 + 551,525 
Ci 8 Vuoretainen, O: The Temperature Field Produced in the Ground by a Heated Slab Laid direct on Ground, 
and the Heat Flow from Slab to Ground. (Acta P. 278/1960) 59 pp. Sw. Kr. 7: 00 
UDC 536.2:697.133:69,025.1 


= 

| 
4 

3 
aa 


q 


